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Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) 
flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a 
circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics 
Research Tunnel (BART). The RANS simulations are compared to a low blowing case 
with a jet momentum coefficient, of 0.047 and a higher blowing case of 0.115. Three 
dimensional simulations of the model and tunnel walls show wall effects on the lift and 
airfoil surface pressures. These wall effects include a 4 % decrease of the midspan sectional 
lift for the C M 0.115 blowing condition. Simulations comparing the performance of the 
Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, 
showing the SST model compares best to the experimental data. A Rotational/Curvature 
Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement 
in the CFD predictions. 


Nomenclature 


BART 

Basic Aerodynamics Research Tunnel 

rh 

Mass flow rate 

CC 

Circulation Control 

NPR 

Nozzle Pressure Ratio, Pt/p< 

CFD 

Computational Fluid Dynamics 

P 

pressure 

GTRI 

Georgia Tech Research Institute 

Pt 

total pressure 

LMP 

Low-Mach number Preconditioning 

Q 

dynamic pressure, \pU 1 2 

NASA 

National Aeronautics and Space 

S 

Reference area, be 


Administration 

U, V, w 

Velocity in x,y,z directions 

RANS 

Reynolds-averaged Navier-Stokes 


respectively 

RCC 

Rotational/Curvature Correction 

V 

Velocity vector, (u,v,w) 


to turbulence model 

X 

Streamwise axis 

RE 

Richardson Extrapolation 

y 

Vertical axis 

a 

angle of attack, deg. 

z 

Spanwise axis 

b 

span (28.0 in.) 



c 

chord length (8.6 in.) 



Cl 

lift force coef., L/(q 00 cb) 

Subscripts 

Cn 

jet momentum coeff. = muj /(qooS) 


Freestream 

c p 

pressure coef., {p — Poo) / (qoocb) 


Jet 

M 

Mach number 

3 
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I. Introduction 


T he design of effective circulation control (CC) configurations require accurate and efficient prediction 
methods. To assess the capability of computational fluid dynamics (CFD) to predict CC airfoil flows, a 
CC workshop was held at NASA Langley in 2004. 1 One of the outcomes from this workshop was that the 
lift performance of a benchmark CC airfoil case was greatly over predicted by almost all of the CFD codes. 
It was also found that the turbulence models played a critical role in the over prediction of lift where the 
jet detachment location on the Coanda surface was typically delayed as compared to experiments. Research 
has shown that by forcing the flow to separate on the Coanda surface, at the same location as observed in 
the experiment, that the CFD will typically match the experimental lift. Therefore, prediction of the correct 
separation location on the Coanda surface is synonymous with a CFD codes ability to predict the correct 
CC airfoil lift. Current research is showing that the RANS turbulence models maybe under predicting the 
eddy viscosity resulting in a CC jet that remains attached longer to the Coanda surface, thus over predicting 
the lift. 2 * 

In an effort to improve CFD codes for CC airfoils, NASA has initiated wind tunnel experiments to 
generate benchmark quality experimental data. This effort focuses on a simple CC airfoil geometry that 
was tested in both the separation control and super-circulation control regimes. Measurements for this 
experiment also focused on characterizing critical boundary conditions required for CFD code validation. 3-5 

Historically code validation comparing 2D CFD to ”2D” wind tunnel test data has been problematic 
because the physics are inherently 3D and typically are influenced by wall effects. While there are simple 
angle of attack and freestream dynamic pressure corrections to approximate wall effects, they are typically 
used for low to moderate lift configurations and are inaccurate for airfoils with very large lift coefficients. To 
remove the ambiguity of comparing 2D CFD to CC wind tunnel data with wall corrections, this study will 
model the tunnel and model geometry in 3D, thereby simulating the actual wind tunnel experiment. The 
results from these 3D simulations will be used to help quantify wall effects and to compare against recent 
benchmark CC experimental data. 4,5 

II. Numerical Modeling Approach 

A. Flow Solver 

The flow field for the circulation controlled airfoil was computed using the flow solver code, OVERFLOW 6, 7 
developed at NASA. This code solves the compressible RANS equations using the diagonal scheme of Pulliam 
and Chaussee. 8 The RANS equations are solved on structured grids using the overset grid framework of 
Steger, Dougherty, and Benek. 9 This overset grid framework allows for the use of structured grids for 
problems that have complex geometries and moving bodies. To improve the convergence and accuracy 
of the flow solution, the OVERFLOW code includes a Low-Mach number Preconditioning (LMP) option. 
The numerical simulations were performed using the parallel version of the OVERFLOW code developed by 
Buning. 10 This code uses the Message-Passing Interface and can run on a tightly-coupled parallel machine or 
a network of workstations. The code distributes zones to individual processors and can split larger individual 
zones across multiple processors using a domain decomposition approach. 

The RANS equations are solved implicitly using the Beam- Warming block tridiagonal scheme with a 3 nd 
order Roe upwind scheme for the inviscid flux terms. The Shear Stress Transport (SST) 11 and Spalart All- 
maras (SA) 12 turbulence models were both used assuming fully turbulent flow including inside the plenums 
and nozzle. A Rotational/Curvature Correction (RCC) model, as implemented in OVERFLOW, was evalu- 
ated for both turbulence models. 13 

B. Grids 

The grids used for the 2D simulations are shown in Fig. 1 and are also used for the cross sections along the 
model span to build up the 3D grids. Figure 2 also shows the tunnel walls and contraction used in the 3D 
simulations. The structured overset grid system for both 2D and 3D were generated using the Chimera Grid 
Tools package. 14 

The airfoil grids have a total of 1241 grid points around the entire surface with 401 of these points 
concentrated on the Coanda trailing edge. There are 101 grid points in the viscous direction away from the 
airfoil surface within a distance of 0.23c, where c is the airfoil chord. The Coanda trailing edge grid has 261 
grid points in the viscous direction, extending 0.23c from the surface. The large number of grid points on 
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Figure 1. Grids used in the 2D simulations and 
for the cross sections along the model span for the 
3D grids. This figure also shows the grids used to 
model the internal plenums. 


Figure 2. Surface grids for the 3D simulation that 
includes the tunnel walls and contraction. The 3D 
simulation also used a symmetry plane about the 
centerline of the tunnel. 


the Coanda surface were used to resolve the evolution of the wall jet and the shear flow region on the outer 
edge of the jet. 

The airfoil body grids are overset onto a single Cartesian background grid that has a uniform grid spacing 
of 0.0058c in all directions and extends 0.47c upstream of the leading edge and downstream of the trailing 
edge. This is the finest Cartesian background grid and is referred to as a level 1 grid. The level 1 background 
grid is then overset onto a group of level 2 background grids that have double the cell size in all directions 
as the level 1 grid. The level of the Cartesian background grids are then doubling in cell size at each level 
until the far field boundaries are approximately 85 chord lengths away from the wing. The 2D simulations 
have approximately 500,000 grid points where 300,000 of these points were used for the body grids. 

The 3D simulations have a symmetry plane at the centerline of the model with viscous grid spacing in 
the span direction at the wall juncture resulting in a total of 142 million grid points. The body grids and 
level 1 background grid, for the 3D case, have 309 points along the model span with 41 points clustered near 
the tunnel wall to resolve the tunnel wall boundary layer and juncture flows. 

III. Wind Tunnel Experiment 

The circulation control experiment for CFD validation was a collaborative effort between Georgia Tech 
Research Institute (GTRI) and NASA Langley. The goal of this collaborative research was to generate 2D 
benchmark quality data of a fundamental circulation control airfoil geometry for CFD validation. 3 4 The CC 
model was developed and tested at GTRI with followup wind tunnel tests at NASA Langley. GTRI per- 
formed extensive testing of the overall performance and parameter characterization, while NASA is pursuing 
additional flow field measurements needed for CFD validation at the NASA Langley BART. 

A. BART Wind Tunnel 

The BART is an open-circuit, atmospheric wind tunnel and was designed for code validation efforts. 15 The 
test section is 28 in. high by 40 in. wide and 10 ft long with a maximum speed of approximately 220 ft/s. 
The BART tunnel is also equipped with wall blowing (on the floor and ceiling), at the test section entrance, 
to decrease the effects of the juncture vortex by energizing the wall boundary layer. 4 Wall pressures were 
also measured on the centerline of the sides walls, opposite the upper and lower surface of the CC model. 

B. CC Model 

A 2D section of the CC model is shown in Fig. 3 and has an elliptic leading edge with a large circular trailing 
edge. This model was chosen for benchmark testing because of its simple geometry and large trailing edge 
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Upper 


Piccolo Secondary 




Figure 3. A 2D section of the CC model showing 
the geometry of the CC-E0020EJ airfoil and the 
upper and lower blowing slots and plenums. 


Figure 4. The 28 in. span CC model installation in 
the NASA BART tunnel mounted vertically with 
pressure taps on the centerline of the side walls. 


radius. The large circular trailing edge is conducive to surface and near-body flow field measurements with 
improved accuracy as compared to CC airfoils with smaller trailing edge radii. The CC model has dual 
blowing slots on the upper and lower portions of the trailing edge; however, in this investigation, the lower 
slot was closed with only the upper slot open. 

The CC model has a chord length of 8.6 in. and was mounted vertically in BART from floor to ceiling, 
as shown in Fig. 4, and has a span of 28 in. The model has an aspect ratio (span/chord) of 3.26 and a tunnel 
height to model chord ratio of 4.6. For the BART tests, the blowing slot was set to a constant height of 
0.020 in (0.00232c). This slot configuration is not optimal for the maximum lift performance but was chosen 
for measurement resolution of the jet at the slot exit. The total pressure of the jet was measured using three 
total pressure probes along the span of the model in the secondary plenum. The total temperature and mass 
flow rate were measured in the air supply line external to the model. 

IV. Results 

A. Slot Jet Profile 

In the CFD simulations, the jet conditions are defined by using an inflow boundary condition at the first 
plenum chamber that sets the total pressure, pt, plenum, and total temperature, T tj plenum- The CFD simula- 
tion then solves the internal flow from the first plenum through a contraction to the second plenum into the 
nozzle and then slot exit. Figure 5 shows streamlines from a 2D simulation revealing the complex flow inside 
the second plenum. At this point it is unknown how much of the internal geometry needs to be modeled. 

To verify the slot velocities from the CFD simulations, detailed hot-wire measurements were made at the 
slot exit with the model on a bench top (i.e. , no external flow) and compared to 2D CFD simulation results. 
A close-up view of the slot exit is shown from a 2D CFD simulation in Fig. 6 with streamlines revealing a 
separation zone behind the slot lip. 

Numerical simulations matching the two blowing conditions of the hot-wire measurements were made and 
compared to the experimental results in Fig. 7 and 8. Initially the CFD simulations included an external flow 
with a freestream Mach number of 0.1. The profiles from the CFD showed a significantly larger jet velocity 
for the same jet NPR and total temperature. Since the bench top measurements did not have an external 
flow, a second CFD simulation using a quiescent flow (Mach=0) condition was performed and resulted in an 
almost exact match of the measured slot velocity profile. Figure 7 shows a low-blowing case at an NPR=1.10 
and Fig. 8 a higher-blowing case at NPR=1.21. The comparisons also show a good match of the slot height 
between the CFD and the experiment. 

Currently the standard definition of the NPR uses the tunnel freestream static pressure, which is measured 
upstream at the tunnel contraction location. These results indicate that this freestream static pressure is 
not characteristic of the effective NPR for the jet velocity and that it should not be used to characterize the 
peak velocity. 

Figures 9 and 10 show the static pressure ratio, p/poo, and total pressure ratio, pt/pt, plenum, at the slot 
exit for the CFD simulation with NPR=1.21. From these figures we see that the total pressure at the core of 
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Figure 5. Streamlines shaded by the velocity mag- 
nitude nondimensionalized by freestream speed of 
sound. These streamlines highlight the complex 
flow inside the second internal plenum. The re- 
sults are from a 2D CFD simulation using an NPR 
of 1.21 with a freestream Mach number of 0.1. 


Figure 6. Contour plot of the velocity magnitude, 
nondimensionalized by freestream speed of sound, 
with streamlines highlighting the flow at the slot 
exit and the wake behind the thickness of the up- 
per slot geometry. This plot is from a 2D CFD 
simulation of the circulation control airfoil with an 
NPR of 1.21 and a freestream Mach number of 0.1. 


the slot is nearly equal to the plenum total pressure with and without external flow. The static pressure at 
the slot exit, with no external flow, is between 0.980 and 0.995 of the freestream value. Figure 9 shows that 
with external flow, the static pressure ratio is greatly reduced at the slot resulting in a higher jet velocity. 
Figures 11 and 12 are contour plots of the static pressure ratio near the slot with and without external flow. 
The airfoil with external flow generates lift accelerating the flow around the Coanda surface, lowering the 
static pressure near the slot exit. Therefore, the increased jet velocity seen in the CFD simulations is a result 
of the external flow reducing the static pressure near the slot, increasing the jet velocity. 

While no slot velocity measurements of the CC model in the wind tunnel were made, it is hypothesized, 
based on the CFD results, that the actual peak jet velocity will be larger than the bench top measurements. 
Based on the results with and without external flow, it is recognized that the jet velocity profiles will also 
need to be measured in the wind tunnel with external flow, in order to validate the CFD at the jet exit. 
Once validated with external flow, the jet profile can be validated on the bench top and only a peak jet 
velocity at the exit will be required in the wind tunnel, assuming the velocity profiles are similar with and 
without external flow. 

Another check on the jet velocity difference with and without external flow is to compare the slot mass 
flow rates from the bench top and wind tunnel experiments for varying NPR values. Figure 13 shows a 
comparison of the mass flow rates between the CFD and experiment, with and without external flow as a 
function of NPR. This figure confirms that both the experiment and CFD have higher mass flow rates at a 
given NPR with an external flow when compared to the bench top cases. The plot also shows that the CFD 
is matching the mass flow rates of the experiment very well for both the quiescent and external flow cases 
and supports the influence of the external flow on the peak slot velocity. 

B. 2D CFD Simulations 

Simulations in 2D were performed in ”free-air” (i.e. , no tunnel walls) with the far field boundaries 85 chord 
lengths away from the body. The simulations matched a freestream Mach number of 0.1, the NPR for the 
internal plenum, and jet total temperature ratio. 

The CFD simulations focused on two blowing conditions from the wind tunnel test, a low blowing case 
at ( 7 ^ = 0.047 and a higher blowing case at C^= 0.115. The low blowing case of ( 7 ^ = 0.047 is considered to 
be in the separation control regime while the higher blowing case of (7^ = 0.115 is at the onset of the super- 
circulation control regime. The airfoil is at 0° for the 2D simulations and no attempt was made to correct 
for the downwash due to the tunnel wall effects. Below we will present the results for the 2D simulations 
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Figure 7. Exit velocity profiles for the circulation 
control model for NPR=1.10. Experimental data 
was measured on the bench top in quiescent flow 
and compared to CFD simulations with and with- 
out external flow. 


Figure 8. Exit velocity profiles for the circulation 
control model for NPR=1.21. Experimental data 
was measured on the bench top in quiescent flow 
and compared to CFD simulations with and with- 
out external flow. 
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Figure 9. Jet static pressure ratio, p/p<x> , profiles 
at the slot exit for the circulation control model for 
NPR=1.21 with and without external flow. 


Figure 10. Jet total pressure ratio, pt /pt, plenums 
profiles at the slot exit for the circulation control 
model for NPR=1.21 with and without external 
flow. 


using the SA and SST turbulence models with and without LMP and RCC. 

1. = 0.047 Case 

A summary of the Cl values for the low blowing case of C M = 0.047 is presented in Table 1. The corresponding 
C p profiles are shown in Fig. 14 for the SA and SST turbulence models with and without RCC and LMP. 
The SA model without RCC predicted a lift coefficient at Cl = 2.03 as compared to the experimentally 
measured value of 1.31 from the model balance. The experiment also had a sectional Cl of 1.35 near the 
3/4 span at j//6 = 0.73 that was integrated from the pressure profile. Note that the sectional pressures near 
the midspan could not be used since this data was missing a critical pressure port on the upper surface of 
the model at x/c = 0.84, resulting in an over prediction of Cl- Therefore, only the sectional Cl near the 
3/4 span location could be used. Adding RCC to the SA model improved the lift results with a Cl = 1.86; 
however, this is still 38% larger than the experimentally measured value. The SST model performed better 
with Cx = 1.58 and was slightly improved by adding RCC where Cl = 1.56. Turning on the LMP had little 
effect on the Cl values and results in a small amount of unsteadiness of the lift convergence history. 

While the suction pressure peak at the trailing edge is close to the experiment, the leading edge suction 
pressure peak is over predicted. It was noted that the 2D simulations are without an angle of attack 
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NPR = 1.21 

Free-Stream Mach = 0.0 



Figure 11. Contour of the static pressure ratio, 
v/Vooi from the CFD simulation for NPR=1.21 
with no external flow, Moo =0. 


NPR = 1.21 

Free-Stream Mach = 0.1 



Figure 12. Contour of the static pressure ratio, 
p/Poo, from the CFD simulation for NPR=1.21 
with external flow, M , oo =0.1. 



Figure 13. Mass flow rate as a function of NPR for the experiment and CFD with and without external flow. 


correction and the experimental data has not been corrected for wall effects. Therefore, any downwash in 
the experiment has not been accounted for in the 2D CFD simulations, which will reduce this leading edge 
suction peak. Simulations modeling the tunnel walls will be discussed in section E below to address the 
effects of the downwash and tunnel walls on the lift coefficient and airfoil pressure distribution. 

2. C M = 0.115 Case 

Results for the higher blowing case of (7^ = 0.115 are shown in Fig. 15 for the 2D airfoil simulation. This 
figure compares the result for the SA and SST turbulence models with and without LMP and RCC. The SA 
model without RCC resulted in a Cl value of 4.55 as compared to the experimental sectional Cl value of 
2.68, which is a 70% over prediction of the lift coefficient. Turning on RCC with the SA model resulted in an 
improved Cl of 3.90; however, this is still 46% larger than the experiment. Using both RCC and LMP with 
the SA turbulence model only reduced Cl by less than 1% and resulted in a small addition of unsteadiness 
to the lift convergence history. Using the SST turbulence model resulted in an improved Cl value of 3.84 
and turning on RCC reduced the value of Cl to 3.34. The SST model with RCC performed best; however, 
the predictions still show a 25% larger lift than the experiment. As in the C IX = 0.047 blowing case, the 2D 
simulations for C^ = 0.115 do not account for downwash on the model or tunnel wall effects by using an 
angle of attack correction. 
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Figure 14. A comparison of the C p profiles using 
different turbulence models from the 2D CFD sim- 
ulations to the experiment for C p = 0.04T with a 
freestream Mach number of 0.1. Experimental C p 
profiles are shown for the mid span of the model 
and at the 3/4 span location. 



Figure 15. A comparison of the C p profiles us- 
ing different turbulence models from the 2D CFD 
simulation to the experiment for C p =0.115 with a 
freestream Mach number of 0.1. Experimental C p 
profiles are shown for the mid span of the model 
and at the 3/4 span location. 


C. 2D Grid Resolution Study 

A grid resolution study was made for the 2D ’free-air’ case for the two C ^ blowing cases using coarse, medium, 
and fine grids. Table 2 summarizes the number of grid points for the three levels of grids where the medium 
grid is representative of the grid resolution of the 2D and 3D simulations. A fine grid was generated with 
a total of 2481 points around the airfoil with 801 points on the Coanda surface. The fine grid also had 801 
points in the viscous direction, normal to the surface that extended a distance of 0.23c. The grid spacing 
at the wall was set to 0.87 x 10~ 7 for the fine grid so that the coarse grid would have y + < 1 at the wall. 
The fine grid case resulted in a total of approximately 1.8 million grid points for the body and background 
grids. The grids for the medium case were generated by removing every other grid point for the body grids 
resulting in a quarter of the grids points as the fine body grids. Similarly, the coarse grids were generated 
by removing every other grid point in the medium body grids. In all three cases the background grids were 
automatically generated by the flow solver where the resolution of the level 1 grid, which encloses the body 
grids, had a grid spacing of 0.0029c for the fine grid case. This background spacing was doubled for the 
medium grid case to 0.0058c and doubled again for the coarse grid case to 0.0116c. 

The grid resolution cases were run using the SST turbulence model with RCC and without LMP. The 
results are given in Table 3 for values of 0.047 and 0.115. The table also shows the Richardson Extrapo- 
lation (RE) results for the medium and coarse grids (RE m-c) and the fine and medium grids (RE f-m) . The 
percentage difference from the RE f-m value is given in the %Cl column. The 0.115 C ^ case shows good 
convergence where there is a 1.26% difference for the medium grid case, which is representative of the grid 
resolution used in the 2D and 3D simulations. The 0.047 C M low blowing case also has good grid convergence 
with the medium grid resulting in a 0.85% difference. 

D. C ^ Calculations 

The CC airfoil performance is characterized as a function of the jet momentum coefficient, C^, and is defined 
as: 


Thrust J ( pjVj ■ h)(Vj • h)dA 

Qoo ^ Qoo ^ 


(1) 


where V) is the jet velocity and n is the unit normal vector at the slot exit. For small scale CC experiments it 
is difficult to measure the jet velocity profile at the slot exit since the slot heights are very small. Therefore, 
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Table 1. Summary of the 2D CFD lift coefficients for the SA and SST turbulence models with and without 
LMP and RC. 





Lift Coef., C L 


LMP 

RCC 

= 0.047 

C „ = 0.115 

Exp. balance 

- 

- 

1.31 

2.58 

Exp. 3/4 span; press. 

- 

- 

1.35 

2.68 


off 

off 

2.03 

4.55 

SA turb. model 

on 

off 

2.05 

4.52 


off 

on 

1.86 

3.90 


on 

on 

1.88 

3.87 


off 

off 

1.58 

3.84 

SST turb. model 

on 

off 

1.58 

3.84 

off 

on 

1.56 

3.36 


on 

on 

1.56 

3.34 


Table 2. Summary of the mesh sizes for a 2D grid resolution study. Body grids consist of a main airfoil, 
circular trailing edge and internal nozzle grids. The total points for the body grids are in the body points 
column and the total number of points for the body and background grids are in the total points column. 


Grid 

Main Airfoil 

Circular TE 

Body Points 

Total Points 

coarse 

421x81 

201x131 

94,317 

168,319 

medium 

841x161 

401x261 

373,492 

507,687 

fine 

1,681x321 

801x521 

1,486,527 

1,795,447 


experimentally the integral in Eq. 1 is difficult to evaluate and in practice the C M values are approximated 
as: 


, = m jUj 

Q°oS 


( 2 ) 


where the thrust is approximated as a product of the mass flow rate and a jet velocity Uj normal to the slot 
exit. The value used for Uj can be a source of much confusion and inconsistency between experiments and 
computations and can have a significant impact on the predicted C values. 

The value for the BART CC experiment was determined by measuring the mass flow rate, to. The 
jet velocity at the slot exit, Uj , was then calculated using the following isentropic relations: 


Uj = 



7 RT t j 



1/2 


( 3 ) 


where ptj is the total pressure measured inside the second plenum and pj is the static pressure at the jet 
exit. This static pressure is typically approximated as the freestream static pressure, p^. However, the 
actual static pressure at the slot exit for CC blowing airfoils is much lower than freestream and results in an 
under prediction of the jet exit velocity. 

Therefore, using the isentropic relations the jet velocity can be calculated by measuring the nozzle pressure 
ratio, NPR=p co /p t j and the total jet temperature. The isentropic relations assume that there are no heat 
transfer and frictional effects. While heat transfer effects are assumed to be small in the experiment, the 
frictional effects (i.e., boundary layers) for a slot of this size are not insignificant. 

A comparison of the C M values between the experiment and CFD are given in Table 4 for the two blowing 
cases. These comparisons show that the C ^ values based on the isentropic Uj, assuming pj ss p^, are 
comparable between the experiment and the CFD using the same approach. The table also shows that the 
mass flow rate per unit span are nearly identical between the CFD and experiment. 
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Table 3. Summary of the lift coefficient results for the grid resolution study. The percentage differences are 
from changes from the Richardson Extrapolation of the fine and medium (RE f-m) grids. 


Grid 

II 

O O 

0.047 

%C L 

II 

3. S) 

o o 

0.115 

%C L 

coarse 

1.5574 

0.99 

3.3198 

2.42 

medium 

1.5596 

0.85 

3.3594 

1.26 

fine 

1.5697 

0.21 

3.3914 

0.31 

RE m-c 

1.5604 

0.80 

3.3726 

0.87 

RE f-m 

1.5730 

- 

3.4021 

- 

Exp. (press.) 

1.35 

2.68 


Table 4. A comparison of C p values and jet exit velocities, Uj , between the experiment and CFD 2D simulations 
for two jet blowing cases. 



CV = 0.047 Case 

C„ =0.115 Case 


Exp. 

CFD 

Exp. 

CFD 

NPR 

1.081 

1.081 

1.208 

1.208 

Tt t0a /Ttj 

0.990 

0.990 

0.987 

0.987 

ih(lbm/ s/ ft) 

0.0437 

0.0435 

0.0693 

0.0680 

C M ( Uj isentropic, Pj & p x ) 

0.047 

0.047 

0.115 

0.114 

C M (integrated) 

- 

0.050 

- 

0.120 

C M (uj mean) 

- 

0.045 

- 

0.109 

C/i ( Uj peak) 

- 

0.055 

- 

0.130 

Uj isentropic, Pj « Poo (ft/s) 

372 

374 

576 

577 

Uj mean (ft/s) 

- 

353 

- 

547 

Uj peak (ft/s) 

- 

436 

- 

659 

Uj peak, isentropic, Pj « pu v (ft/s) 

- 

431 

- 

648 


The actual value of can be computed using 
the CFD solution by integrating Eq. 1 across the jet 
exit. For the low blowing case the actual value 
from the CFD solution is 0.050 as compared to the 
approximated value of 0.047 using the isentropic jet 
velocity. For the higher blowing case the integrated 
from the CFD is 0.120 as compared to the isen- 
tropic value of = 0.114. In both blowing cases, the 
Cfj, values based on measuring the mass flow rate and 
the isentropic Uj resulted in a 5-6% lower value of the 
thrust coefficient as compared to the actual C^. Ta- 
ble 4 also shows that the variation in values base 
on the mean tend to under predict C M while using the 
peak velocity the values are over predicted. Based 
on these comparisons, it is important when compar- 
ing C n values between experiment and CFD results 
that the values of C M are computed using a similar 
approach. Therefore, if the experiment is calculating 
based on a measured mass flow rate and Uj from 
isentropic relations, then the CFD should compute C M 
using the same approach. 



Figure 16. Surface C p near the slot exit for the 
0^=0.115 blowing case. 
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Figure 17. A comparison of the 2D and 3D CFD 
results to the experimental C p profile data for 
C p = 0.047 with a freestream Mach number of 0.1. 
Experimental C p profiles are shown for the mid 
span of the model and at the 3/4 span location. 
The C p profile for the 3D CFD results are from 
the mid span of the model. 


Figure 18. Contour plots of C p on the model and 
tunnel walls with surface restricted streamlines 
for the 3D CFD simulation at C p =0.047. Shown 
is the CC model from the tunnel side wall to the 
mid span of the model where the plane of sym- 
metry boundary condition was imposed. 


Table 4 also shows a calculation of the jet peak velocity using the isentropic relations and the static 
pressure on the back face of the nozzle upper lip, pu p , from the CFD solution. Figure 16 shows a close-up 
view of the CC trailing edge and blowing slot with a plot of C v on the Coanda surface for the CFD and 
experiment. This figure also shows the C p values inside the plenum lower surface and on the external surface 
of the airfoil, upstream of the slot. From this figure it can be seen that while the experimental pressure 
measurements downstream of the nozzle are close to the slot, they are not appropriate for computing the 
peak jet velocity since the pressure increases rapidly downstream of the slot. For the C p = 0.047 CFD 
simulation, the jet peak velocity using pi ip and the isentropic relations is 431 ft/s as compared to the peak 
value of 436 ft/s from the CFD simulation. In this calculation the total pressure in the plenum was used 
assuming no losses in the jet core, where the peak jet velocity is being calculated. For the C p = 0.115 
CFD simulation, the jet peak value using pn p and the isentropic relations is 648 ft/s as compared to the 
peak value of 659 ft/s from the CFD simulation. Based on these calculations of Uj, using the local static 
pressure on the nozzle lip and the isentropic relations, the peak velocity of the jet can be predicted within 
approximately 2% as compared to approximately 13% using the freestream static pressure. Therefore, it is 
recommended for CC experiments to measure the local static pressure when estimating the peak jet velocity 
when comparing to CFD. That being said, it is recognized that this type of static pressure measurement 
near the slot exit can be very difficult to acquire accurately. 

E. 3D CFD Simulations 

The difficulty in comparing 2D ”free-air” CFD simulations to the ”2D” wind tunnel experiments is that the 
flow field for the experiment is inherently three dimensional with wall effects. The wall effects are dependent 
on many factors, one of them is the model size as compared to the size of the tunnel and another is the 
amount of lift generated by the airfoil. For circulation control problems, the goal is to generate very large lift 
coefficients, which consequently increase the wall effects. Another issue is that it is difficult to build models 
with very small blowing slots and to be able to set and measure the slots sizes accurately. Also measuring 
the flow accurately near the slot becomes challenging as the slots decrease in size. This results in a trade off 
between keeping the model to some reasonable size, minimizing wall effects, and making sure the model is 
large enough for accurate measurements at the slot. 

One of the ways to account for the tunnel wall effects is to experimentally estimate an angle of attack 
correction and tunnel q correction. 16 These wall corrections have been shown to work well but become 
difficult to derive for high lift models. Also the validity of angle of attack and q corrections needs to be 
questioned for models with very high lift coefficients and low aspect ratios as is the case for this test data. It 
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Figure 19. A comparison of the 2D and 3D CFD 
results to the experimental C p profile data for 
C/j, = 0.115 with a freestream Mach number of 0.1. 
Experimental C p profiles are shown for the mid 
span of the model and at the 3/4 span location. 
The C p profile for the 3D CFD results are from 
the mid span of the model. 


Figure 20. Contour plots of C p on the model and 
tunnel walls with surface restricted streamlines 
for the 3D CFD simulation at 0^=0.115. Shown 
is the CC model from the tunnel side wall to the 
mid span of the model where the plane of sym- 
metry boundary condition was imposed. 


also becomes difficult to isolate wall effects from modeling issues when CFD and experiments do not match. 

Another way to address the wall effects issue, and the approach taken in this investigation, is to perform 
a 3D CFD simulation of the model and the tunnel walls and actually simulate the wind tunnel experiment. 
While this sounds straight forward, additional challenges like increased cost due to the large number of grid 
points required and the difficulty of modeling the wall juncture vortex accurately arise. 


1 . 


C p = 0.047 Case 


From the 2D simulations, it was found that the SST turbulence model with RCC performed best at matching 
the surface pressures at the jet exit and was therefore used in the 3D simulations. Simulating the model and 
tunnel for C p = 0.047 resulted in a Cl of 1.47 for the entire model with a sectional Cl at midspan of 1.53. 
A summary of the lift results are given in Table 5 showing that the sectional Cl at the midspan for the 3D 
CFD simulation results in a 1% lower value as compared to the 2D CFD. Therefore, the wall effects on the 
midspan Cl are shown to be small for the CFD simulations as might be expected for this low blowing case. 

A comparison of the C p profiles is given in 

Fig. 17 showing a slight difference in the lead- Table 5. Summary of Lift coefficient, Cl, between ex- 
ing edge suction pressure peak between the 2D periment and 2D and 3D CFD simulations. Experiment 

and 3D midspan CFD predictions with no dif- data show " “ from the balance for the entire model and 

a sectional Cl integrated from pressures near the 3/4 

ference at the jet exit and Coanda surface. This span of the model at y / b = 0 . 7 32. 
difference in the leading edge peak C p values is 
attributed to the downwash effect in the 3D sim- 
ulations, reducing the effective angle of attack 
on the model thus producing a smaller leading 
edge suction pressure peak. While the 3D re- 
sults are improved the effects of the walls does 
not improve the over prediction of the lift as seen 
in the 2D CFD cases. Surface C p contours and 
streamlines are shown in Fig. 18 highlighting the 
effects of the juncture flow on the tunnel walls 
and model. 



Lift Coef., Cl 


C p = 0.047 

C M = 0.115 

Exp. balance 

1.31 

2.58 

Exp. 3/4 span; press. 

1.35 

2.68 

2D CFD 

1.56 

3.36 

3D CFD; model 

1.47 

3.09 

3D CFD; midspan 

1.53 

3.21 
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2. C p = 0.115 Case 

The C p = 0.115 case has a much higher Cl and is expected to have larger wall effects as compared to the 
lower blowing case. The 3D CFD simulation for the C p = 0.115 case resulted in a C'l = 3.09 for the whole 
model with a sectional Cl value of 3.21 at the midspan of the model. This case only shows a 4% decrease 
in the midspan Cl as compared to the 2D CFD Cl of 3.36. Therefore, the C p of 0.115 also shows a small 
effect of the walls on the midspan sectional lift coefficient as was seen in the low blowing case. 

A comparison of the C p profiles between the CFD and the experiment is made in Fig. 19. This comparison 
shows a difference in the leading edge suction pressure peak between the 2D and 3D CFD results resulting 
in a lower lift prediction for the latter. This comparison also shows no difference in the C p values at the 
jet exit and Coanda surface for the 3D simulation. However, the 3D CFD is still over predicting the lift 
sectional coefficient at the midspan by 20% showing that the CFD is not completely capturing the correct 
physics on the Coanda surface as in the 2D simulations. 

Figure 20 shows the surface contours and streamlines for the 3D CFD simulation and the effects of the 
extent of the juncture flow on the model and tunnel walls. This figure also reveals some interesting waviness 
of the separation location on the Coanda surface. 


3. Tunnel Wall Pressures 


The tunnel wall pressures were measured on the 
centerline of the walls opposite of the upper and 
lower surface of the CC model. The wall pres- 
sure taps begin near the tunnel contraction and 
extend down the test section and end just before 
the tunnel diffuser. Figure 21 compares the up- 
per and lower wall pressures for the baseline, no 
blowing condition for the CFD and experiment. 

The plot shows very low C p values upstream and 
downstream of the model as expected since there 
is no lift on the model at this condition. There is 
a large jump in the tunnel wall C p values about 4 
chord lengths downstream from the trailing edge, 
which corresponds to a 2 in. wide slot on the ceil- 
ing at x/c= 4 that extends across the ceiling and 
ends 6.5 in. from each of the side walls. This slot 
was used for a 7-hole probe rake to measure the 
downstream wake flow; however, it has resulted in 
contamination of the wall pressures. 

The wall pressures for the two blowing conditions investigated in this paper are shown in Figs. 22 and 23. 
These two figures show how the 3D CFD simulations are matching the entrance C p values of the experiment. 
As a result of the over prediction in the lift, the wall pressure at the same blowing conditions are also over 
predicted. Downstream of the model, near the ceiling slot, it is difficult to compare the wall pressure since 
the CFD is not modeling the open rake slot in the wind tunnel. Figures 24 and 25 show that by subtracting 
the baseline wall pressures data and plotting, A C p = C p —C Pt Baseiine, that the pressure jump associated with 
the open rake slot in the tunnel is removed, improving the comparison with the CFD. Overall the general 
shape of the wall C p profiles are similar between the experiment and the CFD. The peak values do not match 
between the CFD and experiment since the CFD over predicts the model lift. 



Figure 21. Tunnel wall centerline pressures for the 
walls on the upper and lower side of the airfoil with 
no blowing. 


4- Spanuiise Model Pressures 

A spanwise row of pressure taps was placed on the upper surface of the model at x/c = 0.69 in order to 
measure the uniformity of the flow along the span. Figure 26 shows the spanwise pressure profile for the 
experiment and the CFD at two blowing conditions. Since the CFD simulation was performed with a plane 
of symmetry along the centerline of the CC model, the C p distribution has been reflected about the mid 
span of the model making it implicitly symmetric. The base of the CC model is at a span location of z = 0in. 
where the model extends through the wind tunnel floor and is mounted on a force and moment balance. The 
model is mounted on a circular base plate that has an overlapping region on the edge of the plate, reducing 
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x/c (Distance from Model Leading Edge) 


Figure 22. Tunnel wall centerline pressures for 
the walls on the upper and lower side of the airfoil 
for the Cfj, =0.047 blowing case. 


Figure 23. Tunnel wall centerline pressures for 
the walls on the upper and lower side of the airfoil 
for the 0^=0.115 blowing case. 




Figure 24. Tunnel wall centerline A C p values on 
the upper and lower side of the airfoil for the 
=0.047 blowing case. 


Figure 25. Tunnel wall centerline A C p values on 
the upper and lower side of the airfoil for the 
0^=0.115 blowing case. 


the external flow into the test section. The tip of the CC model is at a span location of z = 28 in. and has 
a small gap between the model and the wall. The spanwise profiles from the experiment are non-symmetric 
with an increase in pressure near the base of the model versus the tip. The figure shows that for z <15 in. 
the C p values are higher compared with data at z > 15 in. This increase in C p near the base indicates that 
the model may have a reduction of lift as compared to the other half of the model. The ratio of the sectional 
Cl near the 3/4 span location to the model balance Cl is 1.04 for the experiment at C p of 0.115. The ratio 
of the midspan sectional Cl to model Cl for the CFD is also 1.04, indicating that the loss of lift seen in 
the spanwise pressure for the experiment is very small and is not having a significant impact on the model 
lift. At this time it is unknown as to the source of the non-symmetric Cp distribution and requires further 
investigation. The profiles of C p from near the centerline at z = 15in. to the tip at z = 28in. are similar 
in character between the experiment and the CFD results. The CFD simulations, at matching C p values, 
also show that the CFD significantly over predicts the C p values, which is consistent with the CFD over 
predicting the total model lift. 

The spanwise lift distribution is given in Fig. 27 for the 3D CFD simulations at the two blowing conditions. 
This figure shows that the lift distribution is fairly flat for the two blowing cases over most of the model; 
however, there is a significant loss of lift near the tunnel walls. This loss of lift near the wall results in a 
difference of the lift coefficient at the midspan as compared to the entire model and is more pronounced for 
the higher blowing case, as expected. 
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span location (inches) 

Figure 26. Comparison of the spanwise C p distri- 
bution between the experiment and CFD for two 
blowing conditions. The spanwise pressure ports 
are located at the a?/c = 0.69 station on the upper 
surface of the airfoil. 



Span Location (inches) 

Figure 27. CFD spanwise Cl distribution for C M = 
0.115. 


V. Summary 

R.ANS simulations of a circulation control model used in recent wind tunnel experiments conducted at 
the NASA Langley BART have been performed. The initial CFD investigation started by comparing the 
slot velocity profiles to experimental bench top hot-wire measurements. These comparisons showed good 
agreement to the 2D CFD simulation with a quiescent flow condition at matching NPR values. Using the 
same NPR for a freestream Mach number of 0.1 resulted in a higher peak jet velocity due to the reduced 
static pressure near the exit of the nozzle when lift is generated by the oncoming freestream flow. These 
results showed that the peak jet exit velocity, that would exist in a wind tunnel test or in flight, cannot be 
inferred by direct measurement on the bench top with no external flow and should to be measured in the 
wind tunnel. Based on CFD simulations, the jet profile shape did not change with external flow, suggesting 
that the jet exit profiles can be measured on the bench top and only the peak jet velocity needs to be 
measured in the tunnel, thus reducing the complexity of the wind tunnel test. 

Evaluation of the CFD results showed that the freestream static pressure is not appropriate in calculating 
the jet exit velocity using isentropic relations. It was discovered that a static pressure measurement near 
the slot exit on the upper lip of the nozzle can be used to calculate the peak jet exit velocity within 2%. 
However this measurement can be difficult due to the size of the slot and the rapid change in the static 
pressure downstream of the slot. A static pressure measurement on the back step of the upper nozzle lip 
is experimentally possible and was shown to be a good location for the calculation of the jet peak velocity 
based on the CFD simulations. 

Simulations in 2D were performed to evaluate the SST and SA turbulence models and the effect of RCC 
and LMP on the prediction of the airfoil sectional lift coefficient. The 2D simulations showed that the 
SA turbulence model greatly over predicted the lift as compared to the SST model. In fact, the 2D CFD 
simulation using the SA turbulence model for C M = 0.115 over predicted Cl by 70%. Including the RCC 
to the turbulence models improved the CFD predictions where the SST model with the RCC gave the best 
results. However the SST with RCC still over predicted sectional Cl by 25% without using wall corrections. 
Using LMP did not show a significant effect on Cl but introduced a small amount of unsteadiness to the lift 
convergence history. 

The effects of the tunnel walls were addressed by performing 3D CFD simulations modeling the wind 
tunnel experiment. These simulations showed only a small effect of the tunnel walls on the sectional lift at 
the midspan of the model. The CFD results showed a 4% difference in the midspan sectional Cl values for 
the Cfj, = 0.115 case. Therefore, the wall effects for the two blowing cases investigated in this paper have 
been shown to be very small compared to the over prediction of the lift by the RANS flow solver. 
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